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We study the scissors modes of dipolar boson and fermion gases trapped in a spherically symmet- 
ric potential. We use the harmonic oscillator states to solve the time-dependent Gross-Pitaevskii 
equation for bosons and the time-dependent Hartree-Fock equation for fermions. It is pointed out 
that the scissors modes of bosons and fermions can be of quite different nature. 
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I. INTRODUCTION 

An oscillation of a trapped ultracold gas, the so-called scissors mode, has theoretically been investigated by Guery- 
Odelin and Stringari [J] in connection to superfluidity associated with Bose-Einstein condensation. The mode has 
experimentally been observed in boson gases 0, Q and also in fermion gases Q . The restoring force of the scissors 
mode is originated in the deformation of a trapping potential Trapped ultracold dipolar gases, which have recently 
been attracting experimental and theoretical interests may provide us with opportunities to study a new kind 

of the scissors mode where the dipolar interaction plays an important role as a restoring force. A schematic 
illustration of such a scissors mode is shown in Fig. [1] ^Td\ . Since a tilted configuration (b) is energetically unfavorable 
in a completely polarized dipolar gas, a dipolar restoring force acts to bring the gas to an aligned configuration 
(a) . The aim of this paper is to study the properties of such a scissors mode in a dipolar boson gas and also in a 
dipolar fermion gas. We use the time-dependent mean-field approaches, the time-dependent Gross-Pitaevskii (TDGP) 
equation for a boson gas and the time-dependent Hartree-Fock (TDHF) equation for a fermion gas. We consider gases 
trapped in a spherically symmetric harmonic potential so that pure effects of the dipolar interaction on the scissors 
mode can be elucidated. To solve these equations, we use the harmonic oscillator states with the same frequency as 
the trapping potential, which facilitates the calculation of the matrix elements of the dipolar interaction and allows 
us to treat boson and fermion gases in a similar framework. A disadvantage of our approach is that it is restricted to 
weak interaction regime and to rather small fermion systems. The paper is organized as follows. The formulation is 
given in Sect. II, the results are presented in Sect. Ill and Sect. IV is devoted to summary. 
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FIG. 1: Schematic illustration of a new kind of the scissors mode. The arrows pointing the z axis indicate the polarized 
dipoles. Since a tilted configuration (b) is energetically unfavorable, a dipolar restoring force acts to return the gas to an 
aligned configuration (a). 



II. FORMULATION 



A. Hamiltonian 



We consider trapped dipolar gases of bosons and single-component fermions at zero temperature. We assume 
that the gases are trapped in a spherically symmetric harmonic potential with frequency w and that the dipoles are 
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completely polarized along the z axis. The systems can be described by the hamiltonian 

= ^e^al^aa + i ^ (a/3|w|a'/?')ala]3a^/aa', (1) 

where aj, and are the creation and annihilation operators of a harmonic oscillator state a corresponding to the 
trapping potential, — (^{n + 3/2) with n — Q, 1, 2,...., and {aj3\v\a' P') is the matrix element of the dipolar 
interaction v{r) — {(P /r^){\ — 3z^/r^). Here d is the magnetic dipole moment for a magnetic dipole. For an electric 
dipole, d is given by a{0)8 where a(0) is the atomic polarizability and £ the electric-field strength [ll|. We use the 
units such that h = 1. Assuming that it would be possible to achieve zero scattering length through a Feshbach 
resonance [H, , we neglect the contact interaction described by the s-wave scattering length, as has often been 
done in theoretical investigations of dipolar boson gases [l^,Jl5j. In the case of a single-component fermion gas the 
contact term does not contribute because of Fermi statistics ^ . The harmonic oscillator states with a definite n have 
degeneracy of {n + l)(n-|-2)/2 - fold and are called an oscillator shell [l^. These degenerate levels contain states with 
different orbital angular momentum I and its z component niz. 



B. Time-dependent Gross-Pitaevskii equation 

The TDGP equation for the condensate wavefunction <i>o(''"7^) is given by 

V2 



^ + V^^tir) + v{r~ r>(r', t)dr' 
2m 



Mr,t). (2) 



where m is the mass of a boson, Vcxt(f") — muj'^{x^ + + ^^)/2 and n{r,t) — \^o{r,t)\'^. The wavefunction ^o{r,t) 
is normalized to N, the number of bosons. We expand ^o{r, t) using the time- independent harmonic oscillator states 
as 

$o(r,t)-^C/a(t)0o(r), (3) 

a 

where 4>a{i') is the harmonic oscillator wavefunction. In order to treat boson and fermion gases in a similar framework, 
we derive an equation of motion for a matrix Naa' defined by 

No.o.-{t) = Uo,{t)Uo,,{ty. (4) 

The matrix Naa' satisfies Naa — N. Using Eq. ([2]) and its complex conjugate, we obtain the equation of motion 
for Naa' as [l7t 

iNaa' = ^{ea\Nxa' - Naxexa'), (5) 
A 

where eaa' is given by 

Eaa' = eaSaa' + ^ {aXi\v\a' X2) Nx^X^ ■ (6) 
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Equation ^ can also be derived from the hamiltonian Eq. ([T]) using the time-dependent total wavefunction |<I'(t)) = 
exp[~iHt]\^{t = 0)) and defining Naa' by 

Naa'it)^ imio-laaMt)). (7) 

The time derivative of Naa' is given by 

iNaa' ^ {m\[alaa,H]mt)). (8) 

Using the Hartree approximation for the expectation values of two-body operators on the right-hand side of the above 
equation, 

{^{t)\alal,af,aa\^{t)) ^ Naa'N^p,, (9) 



3 



we obtain Eq. ([5]). Equation ([5]) conserves the total energy Eq of the condensate given by 

Eq = eaNaa 

a 

+ ^ E {af3\v\a'f3')N^^^Nf,,p. (10) 

afla'fl' 

Equation ([SJ also conserves the total number of bosons. The condensate is given as a stationary solution of Eq. ([5]), 
which satisfies 

E(e"A^Aa' - Naxexa') = 0. (11) 
A 

Equation (jlip implies that both e^a' and Naa' can simultaneously be diagonalized. An eigenstate of Naa' which has 
eigenvalue N is the condensate. 

C. Time-dependent Hartree-Fock equation 

The TDHF equation for a fermion gas can be written in a form similar to Eq. ([S]) using the occupation matrix Uaa' 
defined by 

naa'it) = (*(t)|aj,,a„|*(i)), (12) 

where |5'(i)} is the time-dependent total wavefunction. The equation for riaa' is given as 

in^^, = {^malaa„H]\^{t)). (13) 

Using the Hartree-Fock (HF) approximation for the expectation values of two-body operators on the right-hand side 
of the above equation, 

(^'(t)|aj^,a|j,a/3aa|*(t)) « Uaa'TLpp, - nap'Tlpa', (14) 

we obtain the TDHF equation [l^ 

iflaa' = ^^{eaxnxa' " na\e\a>). (15) 
A 



The energy matrix eaa' is given by 



+ ^{aXi\v\a'X2)An\2\^, (16) 

Ai A2 



where the subscript A means that the corresponding matrix is antisymmetrized. The TDHF equation (jl5p satisfies 
the conservation laws of the total energy and the total number of particles [l^. The total energy Ehf in the HF 
approximation is given by 

+ ^ E (a^l«|a'/3')A"-Q'Q»^/3'/3- (17) 

The HF ground state is given as a stationary solution of Eq. (fTSj) . which satisfies 

^Y^{^a\n\a' - na\exa') = 0. (18) 

A 

Equation (IT51) implies that both eaa' and be diagonalized simultaneously. The eigenvalues of either 

1 or 0. 
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D. Scissors mode 



As mentioned above, the condensate of a boson gas and the HF ground state of a fermiongas are given as stationary 
solutions of Eq. ^ and Eq. , respectively. We use the following adiabatic method [13, [3 Ho] to obtain nearly 
stationary solutions of Eqs. (O and (|15l) : Starting from a noninteracting configuration, we solve Eqs. ([5]) and 
using a temporally increasing interaction v{r) x t/T. The noninteracting configuration for a boson gas is the lowest- 
energy harmonic oscillator state and that for a fermion gas consists of the harmonic oscillators states filled up to a 
certain shell determined by the number of fermions. We must use large T so that spurious oscillations coming from 
mixing of excited states are well reduced. We found that the value T = IStt/w is sufficiently large to suppress such 
unphysical contributions. 

In experiments JJ the scissors mode is excited through a sudden rotation of the trapping potential by a small 
ang le. In our approach this is achieved by using the time-dependent one-body operator k X^aa' {'^\QW)<^aO,a' S{t — T) 
[20|, where Q is the angular momentum operator, Q{r) = = —i{yd/dz — zd/dy), and A: is a parameter determining 
the tilted angle and, therefore, the amplitude of the oscillation. The matrix elements of are nonvanishing only 
between the harmonic oscillator states with the same n and I and with the difference in niz by ±1. The initial 
condition for Naa' at t ^ T becomes 

N^^,{T+) - ^(a|e-*^-«|A)iVAA'(7^-)(A'|e*"'=|a'), (19) 

AA' 

where T_ and r_|_ indicate the times infinitesimally before and after T, respectively. Equation (|19l) gives an initial 
configuration of a gas rotated along the x axis. As a consequence, the expectation value of yz is nonvanishing att = T. 
One may think that the operator yz can also be used as the excitation operator of the scissors mode. However, yz is 
essentially the quadrupole operator r'^Y2±i{0, </>) and excites a deformation mode different from the scissors mode. In 
the case of a fermion gas, N\\i in Eq. (IT^ is replaced by n\\i . 

We study the scissors mode in the small amplitude regime and, therefore, expand Eq. (jl9p in powers of k: We take 
terms up to the fourth order of k. The linearized form of the TDGP equation may be referred to as the Bogoliubov 
equations [2l| and that of the TDHF theory is the random-phase approximation (RPA) [l^. In the small amplitude 
regime the physical behavior of collective oscillations are usually investigated using the strength function S{E) |T8l.[22j 
(the dynamic structure factor (23j ) for the excitation operator Q, which describes the transition strength at excitation 
energy E. In our approach, S{E) is calculated as [l^l 

1 

S{E)^— {q{t) -q{T)) sin Et'dt', (20) 
kn Jo 

where q{t) = (L^) and t' = t — T. Since the integration in Eq. (|20|) is performed for a finite interval in our numerical 
calculations, we multiply q{t) — q[T) by a damping factor exp(— ri'/2) to suppress spurious oscillations in S{E). The 
strength function satisfies the following energy-weighted sum rule expressed with the double commutator of H and Q 



ES{E)dE = i($o|[Q,[i/,g]]|$o), (21) 



where |$o) is the ground-state wavefunction. The right-hand side of the above equation for H in Eq. ([T|) and Q = 
is written as d^(<i>o|M^|<i>o)/2, where W is a two-body operator W = 3{{z - z')^ - iv - 2/')^)/!'" - '"'l^- Thus the value 
of the energy-weighted sum for the scissors mode directly depends on the interaction strength in contrast to other 
collective modes such as the breathing mode and the quadrupole mode [l^ [13 ■ For example the right-hand side of 
Eq. ([2T|) for the yz quadrupole mode excited hy Q — yz is given by ($o|2/^ + z'^\^o)/2m. 



III. RESULTS 



A. Dipolar boson gas 

We perform our calculations for a gas of 100 bosons. According to Eq. ([2]), the obtained results appropriately 
normalized by N and w depend only on the dipolar parameter C = Ncf/S^^w, where ^ is the oscillator length 
^ — y/l/muj. The condensate energy Eq calculated using the oscillator states up to the n — 8 shell (solid line) is 
shown in Fig. [5] as a function of C. The results obtained using the oscillator states up to the n = 4 (dotted line) and 
6 (dot-dashed line) are also shown for comparison. More than 93% of the interaction energy {Eq{C = 0) — Eq{C)) is 
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FIG. 2: Condensate energy for A'^ — 100 as a function of the dipolar parameter C — NcP /^^uj. The dotted, dot-dashed and 
sohd hne depict the results obtained using the harmonic oscillator states up to the n = 4, 6 and 8 shells, respectively. 
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FIG. 3: Contour plot of the density n(r) in the x — z plane for N — 100 and C = 5. The white area is the most dense. The 
density has axial symmetry along the z axis and reflection symmetry with respect to the x — y plane. 

estimated to be achieved at C = 5 using the oscillator states up to the n = 8 shell. In the following we present the 
results for C < 5 obtained using the oscillator states up to the n = 8 shell. Our calculations suggest that the dipolar 
gas is stable in the range of the dipolar parameter up until C k, 5, which agrees with other theoretical investigations 
of the ground states and excited states of dipolar boson gases [iM, [l^l- The energy gain shown in Fig. [2] is achieved 
by deformation caused by mixing of higher oscillator states. The condensate has a prolate shape and its deformation 
increases with increasing interaction strength. The density profile at C = 5 is shown in Fig. [31 

The expectation value of is shown in Fig. |4]as a function of time for C = 5. The time evolution slightly differs 
from a harmonic oscillation, indicating a mixture of different frequencies. The strength function obtained from Eq. 
(PU)) (solid line) is shown in Fig. [5] The artificial width used is F/cj = 0.1. The dotted line depicts the unperturbed 
strength distribution obtained from the eigenstates of Naa' (and also e^a' ) at i = r_ . The small difference between 
the two calculations indicates that the condensate-particle correlations such as those included in the Bogoliubov 
equations are small. The small peaks located at higher energies are originated in the mixing of higher oscillators 
states into the condensate. The excitation energy of the lowest-energy peak of the scissors mode is shown in Fig. 
ini as a function of C. The excitation energy is 2ti; at C = and increases almost quadratically with increasing C. 
The fact that the excitation energy is close to 2ui in the small C region means that the trapping potential plays an 
essential role as a restoring force in this region. We found that the excitation energy of the scissors mode coincides 
with that of the quadrupole mode excited by the operator Q — yz (or xz). This may be understood by the fact that 
major components of both modes consist of the transitions from the condensate to the lowest-energy positive-parity 
states with = ±1. The difference is that the scissors mode consists of the transitions with A/ = 0, whereas the 
yz quadrupole mode can have the transitions both with = and 2. The energies of these states with vriz = ±1 
are 2u higher than the energy of the condensate at C = and weakly depend on C, while the condensate energy is 
a rapidly decreasing function of C as shown in Fig. [2] This explains the behavior of the excitation energies of the 
scissors mode (and the quadrupole mode ) shown in Fig. [51 
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FIG. 4: Time evolution of (Lx) calculated using the TDGP equation for C — 5 and k — 0.001. 




FIG. 5: Strength function calculated at C = 5. The dotted line indicates the unperturbed distribution. 

In the following we discuss the energy-weighted sum of the transition strength. The transition strength and its 
energy- weighted sum give useful information to understand the dynamical properties of the condensate '23], though 
they may not have experimental significance because the scissors mode is excited not by an external field but by a 
sudden rotation of the trapping potential. The value of the energy-weighted sum is shown as a function of C in Fig. 
[71 As mentioned above, the sum depends on the interaction strength: It is zero at C = and increases with increasing 
C. The behavior of the transition strength of the scissors mode at small C may be explained by the small mixing of 
the oscillator state with n = 2, I ^ 2 and niz — into the condensate because the scissors mode is mainly excited 




FIG. 6: Excitation energy of the lowest-energy scissors mode as a function of C. 
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FIG. 7: Energy-weighted sum of the transition strength as a function of C. 

through the transitions from the component of the oscihator state with n = 2, I = 2 and ttIz = in the condensate to 
the oscillator states with n = 2, I = 2 and niz = ±1. The small mixing of the oscillator state with n = 2, I = 2 and 
niz = in the condensate gives the vanishing transition strength at small C. On the other hand the yz quadrupole 
mode at small C is mainly excited by the transitions between the harmonic oscillator state with n = 0, Z = and 
rUz = in the condensate and the oscillator states with n = 2, I = 2 and ruz = ±1. Therefore, its energy-weighted 
sum, which is given by ($o|2/^ + z'^l^o) /2m, is almost independent of the interaction strength: The change in the 
energy-weighted sum for the yz quadrupole mode is in fact only 10% when C is increased to 5. The small transition 
strength of the scissors mode implies that the scissors mode may be easily masked by the yz quadrupole mode if there 
is a tiny change in the shape of the trapping potential during the excitation process of the scissors mode. 

B. Dipolar fermion gas 

Since the number of the harmonic oscillator states needed in our approach rapidly increases with increasing number 
of fermions, we can only manage rather small systems. We mainly consider gases with N w 20. Such small N systems 
have often been used for theoretical investigations of dipolar fermion gases ^25.] and may be realized in the array 
of microtraps or optical lattices as discussed in Refs. [25l - [27l |. It is pointed out in Ref. [l^l that in a spherically 
symmetric trapping potential, dipolar fermion gases with the closed-shell configurations where the oscillator states 
up to a certain oscillator shell are fully occupied have small interaction energy and negligible deformation because 
of the cancelation between the attractive and repulsive parts of the dipolar interaction. The numbers of fermions 
corresponding to the closed-shell configurations are 4, 10, 20, 35, 56, and so on. Since such closed-shell gases are 
rather exceptional, we first consider a deformed gas consisting of 18 fermions, which has a simple configuration with 
two fermions removed from the A'^ = 20 closed-shell configuration. Since the dipolar interaction favors a density 
distribution elongated along the polarized axis, a configuration where the oscillators states with the highest jm^j are 
empty has the lowest energy among the configurations with = 18 as long as the interaction is not strong: These 
vacated states are the oscillator states with n = 3, Z = 3 and = ±3. When the dipole interaction is absent, there 
is no restoring force for the scissors mode of the gas because all harmonic oscillator states with the same n and / and 
with different niz are degenerate. Therefore, pure effects of the dipole interaction can be seen in the scissors mode 
of this TV = 18 configuration. The HF energy calculated for such a configuration is shown in Fig. [8] as a function of 
X — (P/u!^^. The dotted, dot-dashed and solid line depict the results obtained using the harmonic oscillator states up 
to the n — 3, 5 and 7 shells, respectively. Since more than 90% of the interaction energy {Ehf{x = 0) ~ Ehf{x = 2)) 
is estimated to be achieved, we present the results obtained using the oscillator states up to the n — b shell. Since 
the noninteracting state is already deformed, mixing of the oscillator states is small in contrast to the boson gas 
considered above. The change in the deformation is also small: The quadrupole moment decreases by 5% when x 
increases to 2. As a consequence, the HF energy almost linearly decreases with increasing %. The density profile 
calculated for A^ = 18 and x = 1 is shown in Fig. [HI 

The time evolution of {Lx) at x = 1 is shown in Fig. [TUl It shows a sinusoidal oscillation pattern. The strength 
distribution obtained from Eq. (j20|) is shown in Fig. 1111 The unperturbed distribution obtained from the eigenstates 
of €aa' at < = T_ is also shown with the dotted line. The difference between the solid line and dotted line is attributed 
to particle - hole correlations. Figure [TT] indicates that the effects of particle - hole correlations are quite small. As 
mentioned above, the transition induced by the operator is allowed only between the oscillator states which have 
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FIG. 8: Hartree-Fock energy as a function of x = /^^uj for N = 18. The dotted, dot-dashed and sohd hne depict the resuhs 
obtained using the harmonic oscillator states up to the n = 3, 5 and 7 shells, respectively. 




FIG. 9: Contour plot of the density n{r) in the x — z plane obtained in HF for A'^ = 18 and x = 1- The white area is the most 
dense. The density has axial symmetry along the z axis and reflection symmetry with respect to the x — y plane. 

the same n and I and differ in rriz by ±1. In the case of — 18 the single-particle transitions from the occupied 
negative-parity states with = ±2 to the unoccupied negative-parity states with to^ = ±3 are dominant and the 
major components of these single-particle states are the harmonic oscillator states with / = 3. 

The excitation energy is shown in Fig 112] as a function of x- The fact that the excitation energy is almost 
proportional to x reflects the linear dependence of the HF energy on x shown in Fig. [H The excitation energy is 
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FIG. 10: Time evolution of (L^) calculated in TDHF for iV = 18, x = 1 and fc = 0.01. 
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FIG. 11: Strength function calculated in TDHF for A'^ = 18 and x = 1- The dotted line depicts the unperturbed distribution. 




FIG. 12: Excitation energy as a function of x for A'^ — 18. 



mainly determined by the difference in the potential energy between the negative-parity states with jm^l = 2 and 
\mz\ = 3. Since the scissors mode is dominantly excited by the allowed transitions between the above-mentioned 
occupied and unoccupied single-particle states, the transition strength weakly depends on x in contrast to the boson 
gas: The increase in the transition strength is 22% when x is increased from 0.1 to 2. Thus the x dependence of the 
energy-weighted sum is almost given by that of the excitation energy in the case of the fermion gas considered here. 

The strength function for the yz quadrupole mode is shown in Fig. 1131 The small peak in the low energy region 
corresponds to the scissors mode. The scissors mode of this fermion gas is energetically well separated from the yz 
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FIG. 13: Strength function for the yz quadrupole mode calculated in TDHF for A'^ = 18 and x ~ ^- 
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FIG. 14: Excitation energy of the scissors mode at x = 1 as a function of A'^. 




FIG. 15: Hartree-Fock energy as a function of x for A'^ — 20. 



quadrupole mode, in marked contrast to that of the boson gas. 

The mass dependence of the excitation energy calculated at x = 1 for the deformed configurations with N = 2, 8, 
18, 33 and 54 is shown in Fig. [M) For a fixed strength of the dipolar interaction the excitation energy of the scissors 
mode gradually increases with increasing N. 

Now we present the results for the scissors mode of a gas with a closed-shell configuration, though it might be 
difficult to experimentally observe the oscillation because of a small deformation of the gas. We consider the TV = 20 
system which has the n — 3 closed-shell configuration. The HF energy is shown in Fig. [TS] as a function of x- Since 




FIG. 16: Time evolution of (L^) calculated in TDHF for iV = 20, x = 1 and fe = 0.01. 
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FIG. 17: Strength function calculated in TDHF for A'^ = 20 and x = 1- The dotted line depicts the unperturbed strength 
distribution. The dot-dashed line indicates the TDHF result at x = 0.5. 



the attractive part and the repulsive part of the dipolar interaction are almost canceled each other, the interaction 
energy is quite small. The deformation of the gas is also small [13]. The HF energy decreases nearly quadratically 
with increasing x because the energy gain is achieved through mixing of higher oscillator states as in the case of the 
boson gas. The time evolution of {Lx) is shown in Fig. [16] for % = 1. The strength function in TDHF (solid line) 
is shown in Fig. 1171 together with the unperturbed distribution (dotted line). The samll difference between the two 
calculations indicates small effects of particle - hole correlations. The scissors mode of the close-shell gas has a broad 
frequency distribution centered around 2a;, which reflects the degeneracy lifting within each oscillator shell due to the 
dipolar interaction. The average of the excitation energies stays around 2uj independently of Xi can be seen from 
Fig. [T7]where the result obtained for x — 0.5 is also shown with the dot-dashed line. The fact that the average of the 
excitation energies is about 2a; implies that the trapping potential plays a role as a restoring force. The scissors mode 
of the closed-shell gas has small transition strengths, as can be understood from the comparison of the heights of the 
strength functions shown in Figs. [TT] and 1171 The energy- weighted sum of the strength function in Fig. [17] is only 
6% of that in Fig. [TT] Fermion gases with the closed-shell configurations have no allowed transitions for the scissors 
mode because all the oscillators states with the same n and / and with different ruz are almost completely occupied. 
The scissors mode can be excited through a tiny mixing of the oscillators states caused by the dipolar interaction. 
This explains the small transition strength of the scissors mode of the gas with = 20 . We also calculated the yz 
quadrupole mode for = 20 and found that it shows an oscillation with frequency 2a; as the usual quadrupole mode 
excited by the operator Y2q{6) [2^. Thus the oscillation pattern of the scissors mode of the closed-shell gas is also 
distinct from that of the yz quadrupole mode. 



IV. SUMMARY 



We studied the scissors modes of dipolar boson and fermion gases trapped in a spherically symmetric harmonic 
potential. We used the harmonic oscillator states to solve the time-dependent Gross-Pitaevskii equation for bosons 
and the time-dependent Hartree-Fock equation for fermions. It was found that the scissors mode of the dipolar 
boson gas shows a nearly harmonic oscillation with quadratically increasing frequency with increasing strength of 
the dipolar interaction. It was pointed out that both the trapping potential and the dipolar interaction act as a 
restoring force of the scissors mode of the boson gas. In the case of the dipolar fermion gas it was clarified that 
the properties of the scissors mode depend on the ground-state configuration. The scissors mode of the gas with a 
deformed configuration has linearly increasing excitation energy with increasing strength of the dipolar interaction, 
reflecting a pure effect of the dipolar interaction as a restoring force. On the other hand the scissors mode of the 
gas with a closed-shell configuration shows a complicated oscillation pattern and has a broad frequency distribution 
centered around twice the frequency of the trapping potential, indicating that the trapping potential also plays a role 
as a restoring force, similarly to the boson gas. It was also pointed out that the scissors mode of the dipolar fermion 
gas can be well distinguished from the yz (or xz) quadrupole mode, in contrast to that of the boson gas which has 
the same excitation energy as the yz (or xz) quadrupole mode. Our study was restricted to fermion gases with simple 
ground-state configurations. It is an interesting problem to investigate the scissors mode of fermion gases consisting 
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of multiple configurations, where the effects of ground-state correlations may have to be considered. 
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